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ABSTRACT 

Aims. Type III radio bursts are intense radio emissions triggered by beams of energetic electrons often associated with solar flares. 
These exciter beams propagate outwards from the Sun along an open magnetic field line in the corona and in the interplanetary (IP) 
medium. 

Methods. We performed a statistical survey of 29 simple and isolated IP type III bursts observed by STEREO/Waves instruments 
between January 2013 and September 2014. We investigated their time-frequency profiles in order to derive the speed and acceleration 
of exciter electron beams. 

Results. We show these beams noticeably decelerate in the IP medium. Obtained speeds range from ~ 0.02c up to ~ 0.35c depending 
on initial assumptions. It corresponds to electron energies between tens of eV and hundreds of keV, and in order to explain the 
characteristic energies or speeds of type III electrons (~ 0.1c) observed simultaneously with Langmuir waves at 1 au, the emission of 
type III bursts near the peak should be predominately at double plasma frequency. Derived properties of electron beams can be used 
as input parameters for computer simulations of interactions between the beam and the plasma in the IP medium. 

Key words. Sun: radio radiation - Sun: particle emission - plasmas - methods: data analysis 


1. Introduction 


Type III radio bursts are produced by beams of suprathermal 
electrons accelerated near the Sun’s surface during solar flares 
( Wild|1950| . These exciter beams propagate outwards from the 
Sun along an open magnetic field line in the corona and the in¬ 
terplanetary (IP) medium at large distances beyond 1 au, where 
suprathermal electrons can be detected in situ by spacecraft. 
They locally produce a bump-on-tail instability that leads to a 
generation of electrostatic Langmuir waves near the local plasma 
frequency / pe . The conversion of Langmuir waves into escap¬ 
ing electromagnetic emission is believed to be due to nonlin¬ 
ear interactions involving electromagnetic waves: type III radio 
bursts near / pe (the fundamental, E-component) or near 2/ pe (the 
harmonic, //-component), which can be observed remotely ( Lin j 
1970, Gumett & Anderson|1976) l. This conversion involves the 
wave-wave interactions of Langmuir, ion acoustic and electro¬ 
magnetic waves ( |Ginzburg & Zhelezniakov||1958[ ). Since elec¬ 
tron beams propagate away from the Sun, the characteristic 
emission frequency decreases to consecutively lower frequen¬ 
cies corresponding to a decrease in / pe . 


Although we can distinguish the E- and //-components from 
metric to decametric wavelengths ( i.e. of coronal origin), it is 
difficult to determine the observed component of type III bursts 
at longer wavelengths, which are generated in the IP medium. 
However, we can recognize the F- and //-components for cases 
where type III triggering electron beams intersect the spacecraft, 
those associated with Langmuir waves observed in situ ( [Kellogg] 
|1980)>. Nonetheless, such events occur rather seldom. Although 


IP type III bursts are generally assumed to be emitted at the 
harmonic of the electron plasma frequency (Reiner et al. 1992 


|Mann et al.|1999| >, we will discuss presence of both components 
separately. 

Since early observations, the characteristic exciter speeds as¬ 
sociated with type III bursts have normally been found in the 
range 0.1 - 0.6c. [Dulk et aT| ( j 1987} ) investigated both onset 
and peak times of type III bursts associated with measurements 
of Langmuir waves and suprathermal electrons by the ISEE-3 
spacecraft assuming presence of the E-component alone ( |Scarf| 
|et al.|1978)|Anderson et al.|1978) . They showed the exciter beam 
speed of 0.14c and 0.07c using onset and peak times, respec¬ 
tively. Hoang et al. ( jl994| ) have studied onset times of the type 
III bursts observed by the URAP instrument onboard the Ulysses 
spacecraft ( Stone et al.|1 992) using the same approach as Dulk 
|et al.| ( [T987 l. |Hoang et al.| ( |1994) concluded that the exciter beam 
speed ranges from 0.04c to 0.13c. However, the variation in these 
speeds proved to be more challenging. Fainberg & Stone ( 1970) 1 
have suggested that the beams are propagating over 1 au with 
little deceleration. Later, Poquerusse et ah] ( | 1996} ) suggest that 
deceleration occurs before about 0.1 au, and that there is no 
significant correlation between the initial and final speeds. The 
STEREO provides radio data with significantly better resolution 
than used in the above-mentioned studies. 


In this paper, we present statistical results on speed and ac¬ 
celeration of exciter beams of IP type III radio bursts observed 
by the Solar TErrestrial RElations Observatory (STEREO). In 
Section I 2 TI we focus on the instrumentation. In Section [2721 we 
describe our data set. Then we show a detailed analysis of one 
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IP type III radio burst to illustrate the methods we used (Sec¬ 
tion 2.31. We present statistical results for exciter speeds and ac¬ 
celeration in Section l2~4l 


2. Observations and data analysis 

2.1. The STEREO spacecraft 

STEREO consists of two identical spacecraft orbiting the Sun 
and providing a unique stereoscopic view of the solar-terrestrial 
system ( |Kaiser et al.||2008} . The STEREO/Waves instrument 
records electric field fluctuations by three monopole anten¬ 
nas ( |Bougeret et al. 2008 Bale et al. 2008} . For analysis of 
type III radio bursts we use data recorded by the Low Fre¬ 
quency Receiver (LFR, 10 kHz - 160 kHz in 32 logarithmically 
spaced frequency channels) and High Frequency Receiver (HFR, 
125 kHz - 16.025 MHz in 319 linearly spaced frequency chan¬ 
nels). The time resolution of both receivers is 38 seconds. With 
LFR and HFR we can investigate radio emissions that originate 
somewhere from 2 Rs up to the 1 au of the Sun when assuming 
an average electron density model in the solar wind (Cecconi 
et a l,|2008| l. However, we are unable to accurately calculate fre¬ 
quency drifts of type III radio bursts above 1 MHz owing to low 
time resolution of the STEREO/Whves instrument. Therefore we 
limit our analysis to below 1 MHz. Both receivers can operate 
either in the three-monopole mode or in the dipole-monopole 
mode. We analysed the data only in the three-monopole mode. 
The LFR and HFR have been operating in this mode since Jan¬ 
uary 2013 and May 2007, respectively. We subtracted receiver 
background levels from the data to improve the signal-to-noise 
ratio. These levels were calculated as median values over a long 
period of 15 days for each channel-antenna configuration sepa¬ 
rately. More details on method used can be found in|Krupar et al.| 
|20lQl[20l2l[20T4bl l. 


2.2. Data set 

We performed a statistical analysis of 29 type III radio bursts ob¬ 
served by STEREO/Waves between January 2013 and Septem¬ 
ber 2014. Only intense, simple, and isolated cases (i.e. without 
overlapping) with a low frequency cutoff below 100 kHz have 
been included in our data set. These constraints resulted in rel¬ 
atively few suitable events for time - frequency investigation. 
However, this helped to select radio bursts a with smooth mono¬ 
tone relation between peak time and frequency. 

2.3. September 27, 2013 type III radio burst 

We show an analysis of an IP type III radio burst on Septem¬ 
ber 27, 2013 to illustrate the methods used. Figures 
display the flux density S measured by LFR and HFR onboard 
STEREO-A and STEREO-B. Both spacecraft detected the sim¬ 
ple and isolated type III radio burst with a starting time of about 
14:20 UT. The low frequency cutoff is 60 kHz and 80 kHz 
for STEREO-A and STEREO-B, respectively. We excluded fre¬ 
quencies above 1 MHz because of the insufficient time resolution 
of the receiver. During this event STEREO-A was 140° west of 
a Sun-Earth line at 0.97 au from the Sun, whereas STEREO-B 
was 146° east and 1.05 au of the Sun. 

We investigated times corresponding to peak fluxes since on¬ 
set times are difficult to accurately identify owing to density 
variations that result in fluctuations of the quasi-thermal noise 
(Meyer-Vernet & Perche fl989} . Frequency channels have been 
converted to radial distances from the Sun considering both the 
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Fig. 1 . Analysis of measurements recorded from 14:15 to 16:00 UT on 
September 27, 2013: the flux density 5 for STEREO-A (panels a, and 
b) and STEREO-B (panels c, and d). 


F-component (ty) and the //-component ( r h) using the elec¬ 
tron density model in the IP medium by Sittler & Guhathakurta) 
( 1999} . It allows the obtained times to be corrected for a light 
travel time frequency by frequency from the source to the space¬ 
craft with a radial distance from the Sun (rs/c): 


tAf) = t(f)- 


rs/c ~ r,(f) 


(1) 


where / represents frequency, c is a speed of light, and i stands 
for the F- and //-components. This correction assumes that radio 
sources are located on the Sun-to-spacecraft line. The validity of 
this assumption has been tested using computer simulations in 
Appendix A. 

Previous studies of exciter beam speed were dedicated to 
events associated Langmuir waves observed in situ (]Dulk et al.j 
|1987| |Hoang et al.|1994| |Robinson & Caims|1994| l resulting in 
different assumptions on radio source locations. Plasma param¬ 
eters measured onboard were used to estimate the Parker spi¬ 
ral geometry together with the density profile along it ({Parker] 
1958} . However, this approach does not consider a propagation 
of plasma parameters between the Sun, so the spacecraft located 
at 1 au. It can take up to several days for the solar wind to 
travel between the corona and the spacecraft and both plasma 
speed and density, used for the Parker spiral reconstruction, may 
differ significantly from measurements in situ. Because of this 
drawback and since our events were not associated with Lang¬ 
muir waves, we assume that the radio sources lies on the Sun- 
to-spacecraft line, so the average electron density model in the 
ecliptic was used. 

A comparison between corrected times and frequencies con¬ 
verted into radial distances allows speed and acceleration of ex¬ 
citer beams to be derived. However, this simple method does 
not take several important ambiguities into account. The first is 
that we are unable to distinguish the F- and //-components of IP 
type III radio bursts. The distance re is typically twice greater 
than rp resulting in about twice greater beam speed and acceler¬ 
ation for the //-component than for the F-component. Therefore 
we discuss our results considering the presence of both compo¬ 
nents throughout the paper. Moreover, triangulated radio sources 
of IP type III radio bursts lie at considerably (five times for the 
F-component and three times for the //-component) larger radial 
distances from the Sun than the electron density models predict 
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Fig. 2. Analysis of measurements recorded from 14:15 to 16:00 UT 
on September 27, 2013: radial distances from the Sun vs. corrected 
times for STEREO-A (in red) and STEREO-B (in blue) assuming the 
F- (panel a) and //-component (panel b). Dotted lines represent results 
of the polynomial model (equation[3]>. Dashed lines show the power law 
model (equation[6]>. Obtained parameters of both models are shown on 
the top. 


(Steinberg et al. 11984 , Krupar et al. |2014a i. It indicates that we 
observed only scattered images of real sources that are appar¬ 
ently very extended. We thus conclude that the presented method 
provides us with a lower bound estimate of the exciter speed. 

Figures |2}t and [2jr display radial distances from the Sun (r F 
and r H ) as a function of corrected times (/ F and f H ). We used two 
models to fit the r x vs. t\ dependence. Initially, we considered a 
constant acceleration, when the speed (vi) is a linear function of 
time (/): 


>'i0i) - voi + a,?i, (2) 

where voi is a beam speed at t, = 0 (an initial speed in [kms -1 ]), 
and ci\ represents the beam acceleration in [kms 2 ], The corre¬ 
sponding relation between r, and /, is then 


n(h) = 


DO + t'Oih + 



(3) 


where r m is an initial distance from the Sun at t, = 0. We de¬ 
note it as "the polynomial model" which represents the simplest 
approach for estimation of both initial speed and acceleration. 

The second model assumes that a beam speed v; has a power 
law dependence on the radial distance r, from the Sun as 

v,(D) = aifo/RsA (4) 


where a x represents the beam speed at 1 Rs (a power law initial 
speed) and /J; denotes power-law index, so that positive values of 


Pi represent acceleration. Since we measure r; vs. t, we need to 
solve the first-order nonlinear differential equation: 


dr x 

dt\ 


ffiln/RsA 


(5) 


By separating variables we obtain the following relation between 
r, and t\. 

n(fi) = [(1 -Aifeh+ C,)R S ^]^, (6) 


where C\ is the constant of integration. Since both numerical 
simulations ( |Kontar|2001 1 and theoretical relationships suggest 
that the velocity should vary as a power of the distance of the 
source from the Sun ( jRobinson||1992| , we denote it as "the power 
law model". The acceleration a\ for this model can be obtained 
by differentiation of equation (4): 

0 , (n)= * =A jAr-L* ,7, 

dt Pl '\r s / R s dt R s \R S / 


We performed a polynomial regression of the second or¬ 
der to estimate type III exciter beam parameters according to 
the polynomial model (equation |3J. Calculated initial speeds 
are greater at STEREO-A than at STEREO-B. This differ¬ 
ence is probably related to the time correction method. The 
stronger signal at STEREO-A suggests that the electron beam 
propagates more towards STEREO-A than STEREO-B. The 
time correction method thus underestimates the calculated speed 
slightly at the latter (Appendix A). The E-component assump¬ 
tion (panel a) provides us with slower speeds when compared 
to the //-component as it can be expected. Our results indi¬ 
cate that the electron beam decelerates (a F ~ -7 km s 2 and 
«h ~ -12 km s ' 2 ) depending on an assumption on the observed 
component. For the power law model, we used a nonlinear least 
squares fit of the spacecraft data with equation (|5J. The param¬ 
eters of the model obtained with a % 2 goodness-of-fit statistic 
weighted by the measurement error are shown in Figure [2] Neg¬ 
ative values of p confirm that the electron beam decelerates. 

Both models fit the data very well. We observe a discrepancy 
between the two models for larger distances from the Sun when 
the polynomial model predicts greater deceleration when com¬ 
pared to the power model. Since exciter beams typically travel 
far beyond 1 au, the exciter power law model seems to be more 
applicable for larger distances from the Sun. However, exciter 
speeds farther from the Sun, and the initial speeds (v;o, a;), where 
we do not have any data, are dubious since they represent an ex¬ 
trapolation of both models. 


2.4. Statistical results 

We performed the above described propagation analysis on 29 IP 
type III radio bursts observed by STEREO /Wuves between Jan¬ 
uary 2013 and September 2014 case by case (Section 2.2). 

Figures [3jr - [3jl display histograms of the acceleration co¬ 
efficient A; and p x . Median values of the acceleration for the F- 

component a F are-4 km s 2 while the acceleration for the 

//-component a H is about 50 % larger than a F . However, both 
Py and /in are comparable for the power law model (~ -0.35). 
Figures [3£ & [3]f show the relation between the acceleration and 
the radial distance from the Sun for both models using median 
values from panels above. While the acceleration is obviously 
constant for the polynomial model (Figure [3^), it decreases for 
the power law model with a radial distance (Figure [3} ). The ac¬ 
celeration for both models are comparable between 5 Rs and 
20 R s . 
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Fig. 3. Panels a - d: Histograms of the acceleration coefficient a (the polynomial model, equation^ and the acceleration coefficient (the power 
law model, equation[6| assuming the F- and //-component for STEREO-A (in red) and STEREO-B (in blue). Dashed lines indicate median values. 
Black dotted lines show zero values. Panels e & f: Acceleration vs. distances for the polynomial and power-law models for STEREO-A (in red) 
and STEREO-B (in blue). Dotted and dashed lines represent F- and //-components, respectively. 


Figure [4] shows results on exciter beam speeds v, vs. radial 
distance r, for all events. Curves represent only the data that have 
been measured by STEREO (i.e. above the low frequency cutoff 
case by case. Figures |4jt - |4ji). Calculated speeds range from 
~ 0.02c to ~ 0.22c and from ~ 0.03c to ~ 0.35c for the F- and 
//-components, respectively. From individual kinematic curves, 
we calculated median beam speeds with the 25th and the 75th 
percentiles (Figures [4j -|4ji). The statistical curves we obtained 
contain jumps due to the termination of some bursts. Median val¬ 
ues of the speeds decrease from ~ 0.09c to ~ 0.04c and from 
~ 0.16c to ~ 0.09c for the F- and //-component, respectively. It 
suggests that the exciter beams of IP type III radio bursts decel¬ 
erate in the solar wind. However, other interpretations of this sta¬ 
tistical results can be proposed: the observed deceleration could 
be an apparent observational effect caused by stronger scatter¬ 
ing and thus larger refractive index of radio waves at lower fre¬ 
quencies. The important role of scattering of IP type III bursts 
has been demonstrated via the Monte Carlo simulations by |The-| 
[jappa et alT]( |2007j >. The power law model provides us with more 
reasonable beam speeds when extrapolated to 1 au, while the 
polynomial model predicts a beam halt before 1 au which is not 
generally observed. 


3. Conclusions 

We investigated the properties of 29 simple and isolated type 
III radio bursts observed by the two STEREO spacecraft. The 
F- and //-components have been considered using two different 
models for estimating the exciter beam speed and acceleration. 
We show a detailed analysis of one event from our data set as an 
example. The presented method provides us with a lower bound 
estimate of the exciter speed. 

Our results suggest that exciter beams decelerate in the IP 
medium (Figure |5J. The collisional losses cannot account for 
this effect: the collisional mean free path of electron with kinetic 
energy ~ 10 keV in typical solar wind plasma is higher than 
1 au and the cumulative effect of collisions on the beam energy 
is negligible. The deceleration can be explained by the energy 
loss of electron beams due to instabilities producing Langmuir 
waves and the effect of solar-wind inhomogeneities on these 
waves. This loss of the energy is a two-step process. First, Lang¬ 
muir waves are generated by the front the beam. Second, these 
wave are refracted in decreasing plasma density towards high 
wave numbers (lower phase velocities), thereby reducing the 
level of waves that can be re-absorbed by the beam (e.g. Kontar 
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Fig. 4. Panels a, b, c, and d: Exciter beam speeds v, vs. radial distance r ; for all events considering the polynomial model (equation [ 3 } and the 
power law model (equation[6) for STEREO-A (in red) and STEREO-B (in blue). Panels e, f, g, and h: Median values (solid lines) and 25th/75th 
percentiles (dotted lines) of data in panels a, b, c, and d. 


|2Q01) >. Large scale simulations of electron transport show that 
a large part of the initial energetic electron energy is lost via 
Langm uir waves due to plasma inhomogeneity (Reid & Kontar 


■4 km s - to 


2013). The obtained acceleration a, ranges from - 
~ -7 km s 2 . Derived exciter speeds vary between ~ 0.02c and 
~ 0.35c, which corresponds to electron energies from tens of eV 
up to hundreds of keV (Figure[4), where the higher speeds corre¬ 
spond to the //-component. A comparison with energies of elec¬ 
tron beams detected in situ suggest that the //-component pre¬ 
dominantly occur in our data set ( ILin et al.[1996) >. However, the 
deceleration could also be related to stronger scattering at lower 
frequencies. Presented properties of type III exciter beams can 
be used as input parameters for computer simulations of beam - 
plasma interactions in the IP medium. 


Appendix A: Time correction method validation 

We validated our time correction method via simulations of elec¬ 
tron beams propagating away from the Sun (Figure [AT) . Simu¬ 
lated beams travel with a given constant speed along the X axis 
while emitting radio emissions at frequencies between 40 kHz 
and 1 MHz (i.e. at a frequency range of STEREO/Waves used in 
this study). We have considered the F- and //-components using 
the electron density model by |Sittler & Guhathakurta| ( |1999) . A 
modeled spacecraft lies on a circle around the Sun with a radius 
of 1 au in the ecliptic (Z=0), so its position can be parametrized 
by angle (j) = arctan E. We applied the time correction method 
(Eq. 0 to observed times by the modeled spacecraft. Our results 
show that differences between initial and derived beam speeds 
are negligible for small angles. For a mean separation angle be- 
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Fig. A.l. Simulated beam speeds vs. deviation angle A0. Colors denote 
beam speeds (red: 0.2c, green: 0.14c, and blue: 0.07c). The fundamental 
and harmonic components are represented by solid and dashed lines, 
respectively. Horizontal lines indicate initial speeds. Dash-dotted lines 
show maximal and mean separation angles between the two STEREO 
in our data set. 

tween the two STEREO spacecraft in our data set (<f> = 50°), this 
difference ranges from 1 % to 7 %. The difference at <p - 100° 
varies between 7 % and 20 %. Moreover, variations between the 
F- and //-components are minor. 

We also investigated the role of Parker-spiral propagation on 
estimates of speed and acceleration (Figure |A.2[ ). We analyzed 
the data from STEREO-A from September 27, 2013 assuming 
the /’-component (Section |2.3[ ). The emitting energetic electrons 
were assumed to travel along four Parker-spiral lines with re¬ 
spect to spacecraft position. The parameters obtained for various 
time corrections are not significantly different. For a comparison, 
results with no time correction are shown in black. We thus con¬ 
clude that the time correction method is applicable to our data 
set and that it provides us with reasonable results for electron 
beam dynamics with typical speeds in the solar wind. 
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Fig. A.2. Analysis of measurements recorded by STEREO-A from 14:15 to 16:00 UT on September 27, 2013 assuming the F-component: Panel a): 
radial distances from the Sun vs. uncorrected (in black) and corrected times assuming Parker-spiral propagation (in blue, red, purple, and orange). 
Panel b): Parker-spiral paths and the spacecraft location (in green). 
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